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In QCD one can change the representation of the gauge group for quarks and/or the gauge group itself. 
Examples of such generalizations are: (a) supersymmetric Yang-Mills theory with gauge group SU(2) or SU(3); 
(b) QCD with SU(2) colour and quarks in the adjoint representation; (c) QCD with three equal-mass quarks. 
Experience and perspectives in simulations of these theories by applying the two-step multi-bosonic algorithm 
is described. As an example of an interesting physical problem in these theories the expected spontaneous CP- 
violation in QCD with three equal-mass quarks is discussed. 



1. INTRODUCTION 

Quantum chromodynamics is most probably 
the correct theory of strong interactions. It be- 
longs to a class of relativistic quantum field theo- 
ries in four dimensions which are believed to de- 
scribe a wide range of physical phenomena con- 
sistently from zero to infinite energies. 

Examples of generalizations of QCD are given 
in the abstracts. Numerical simulations of these 
theories can give important non-perturbative in- 
sight in low energy dynamical features as confine- 
ment and chiral symmetry breaking. Restricting 
numerical studies to the quenched approximation 
is obviously unsatisfactory. In QCD, besides dy- 
namical simulations with two equal mass quarks, 
the investigation of 2 1 light dynamical quark 
flavours is very important because it is a good ap- 
proximation of the physically relevant u-,d- and 
s-quark system. 

A fermion simulation algorithm which is able 
to deal with QCD-like theories for small fermion 
masses in sufficiently large volumes is the two- 
step muhi-bosonic (TSMB) algorithm In 
a recent large-scale simulation ^-^J this algo- 
rithm has been used for the numerical study of su- 
persymmetric Yang-Mills theories which contain 
"gauginos" i.e. Majorana fermions in the adjoint 
representation of the gauge group. The algorithm 
could cope with a potentially difficult situation 
where massless gauginos are playing an impor- 
tant dynamical role in reasonably large physical 
volumes. It is an interesting question whether 



this experience can be extended, for instance, to 
QCD with 2 1 light dynamical quarks. 

2. FERMION ALGORITHM 

The TSMB algorithm is based on the multi- 
bosonic representation of the fermion determi- 
nant 01 : 
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Here Q is the fermion matrix and Nf degenerate 
flavours are considered. The polynomial Pn{x) 
is an approximation of the function x~°', {a = 
Nf/2) in an interval [e. A] covering the spectrum 

of gtg: 



lim Pn{x) 



for X G [e, Al 



(2) 



The inverse determinant of P„ in ([l]) can be 
represented by n bosonic "pseudofermion" fields. 
Since for a good approximation a large number 
(n = 0(1000)) of pseudofermion fields would be 
necessary the direct application of eq. (|^) to nu- 
merical simulations is not practical. In TSMB a 
two-step approximation by a product of two poly- 
nomials is introduced 
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The multi-bosonic representation is applied for 
the first polynomial Pr}^ and the determinant 
of the second one is taken into account in a 
stochastic "noisy correction" step. For very small 
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fermion masses it is advantageous to perform a 
final "reweighting correction" in the process of 
evaluating the expectation values. (For a com- 
prehensive summary of algorithmic details and for 
references see Q.) 

The DESY-Miinster collaboration collected ex- 
tensive experience on the behaviour of the TSMB 
algorithm in an application to supersymmet- 
ric Yang-Mills (SYM) theory [|-|. The total 
amount of computer time used in last year cor- 
responds to ^ 0.01 Tflopyear (sustained perfor- 
mance with 64 bit precision). The largest lattice 
was 12^ • 24. In QCD units obtained from the 
string tension the spatial extension of this lattice 
is ~ 1 fm. The statistics was relatively high: sev- 
eral thousands independent gauge configurations 
per simulation point. Since a Majorana fermion is 
equivalent to half of a Dirac fermion, the gaugino 
corresponds to 

Nf/2^1/4: . (4) 

An important factor in the extrapolation to 
QCD is the dependence on the number of flavours 
Nf . It turned out that both execution time of a 
sweep and autocorrelation of important physical 
quantities are roughly proportional to the num- 
ber of pseudofermion fields ni. This has to be 
chosen such that the acceptance rate in the noisy 
correction step stays high enough. (In our simu- 
lations we always had 85-90% acceptance.) This 
can be achieved if the deviation norm satisfies 

< 5 • 10-^ . (5) 

On the basis of this it is possible to estimate the 
dependence on the number of flavours because, in 
the limit of a small lower limit of the approxima- 
tion interval e ^ A, the relation 

^„(a;e^O,A=l)= " ^ ^ (6) 
n + 1 + a In 

gives a good estimate. This shows that the num- 
ber of necessary pseudofermion fields n\ is ap- 
proximately proportional to N ^ and the cost for 
obtaining a new independent gauge configuration 
to N). 

Let us note in this context that in a future QCD 
simulation one would like to have, for instance, a 
lattice spacing a ~ (4GeV)~^ and a light quark 



mass rriq ~ 4MeV implying avriq ~ 10""^. This 
would correspond to e ~ 10~^, a value which has 
also been exploited in SYM simulations Q . 

Most of the experience on the TSMB algo- 
rithm is concentrated up to now on Wilson 
fermions. The implementation of "moderately 
improved" actions, as e.g. the clover action of 
Sheikholeslami-Wohlert, seems possible with rea- 
sonable effort. Recently the algorithm has also 
been applied for staggered fermions 0. 

An important question for numerical simula- 
tions of odd numbers of quark flavours with Wil- 
son fermions is the effect of the sign of the fermion 
determinant. In eq. (|l|) the absolute value of the 
(real) determinant appears therefore the sign has 
to be taken into account as a reweighting factor 
during the evaluation of expectation values. This 
is a potential problem in at least two respects: the 
computation of the determinant is only possible 
on very small lattices and even if the sign of the 
determinant is known the reweighting with it can 
in principle cause a serious loss of accuracy due to 
cancellations. In SYM the problem is similar but 
the role of the determinant is played by the Pfaf- 
fian resulting from the Grassmann integration for 
Majorana fermions. 

Concerning the first problem a satisfactory so- 
lution is to monitor the sign of the determi- 
nant (or Pfaffian) by the spectral flow method 
The determinant (or Pfaffian) has a positive 
sign at small enough hopping parameters (say, 
K < 0.125) and one can count the number of sign 
changes up to the required (dynamical) hopping 
parameter by monitoring the flow of the eigen- 
values of the hermitean fermion matrix Q = 75Q 
across zero. The problem of strong cancellations 
among contributions with opposite sign remains 
a potential danger. 

In SYM with SU(2) gauge group it turned out, 
however, that the configurations with positive 
Pfaffian strongly dominate Q . In fact for smaller 
values of the hopping parameter K no configu- 
rations with negative PfafRan were found. Neg- 
ative signs at the percent level occurred only for 
K-values beyond K = Kq corresponding to zero 
gaugino mass. 

The sign of the determinant in QCD with 2 ® 1 
light flavours was up to now not investigated. 
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However, the s-quark is presumably not light 
enough for this problem to become serious. 

3. THREE-FLAVOUR QCD 

QCD with three equal mass quarks is not re- 
alized in Nature but from the theoretical point 
of view is rather interesting. An important fea- 
ture is the possibility of spontaneous CP-violation 
at negative quark mass predicted long ago by 
Dashen ^ and Witten |^ . This is a very interest- 
ing theoretical issue especially in view of recent 
and near-future CP-violation experiments in B- 
factories. In the lattice literature the question of 
spontaneous CP- and flavour-symmetry violation 
has been extensively studied but mainly in case 
of A/ = 2 flavours where it can also be 

just a lattice artifact [|12|. A similar phenomenon 
also occurs in other models, for instance, in SYM 
theory with SU(3) gauge group (sec |^). 

Spontaneous CP-violation in QCD with three 
flavours is predicted by the analysis of the low 
energy effective chiral Lagrangian 



8 dependence in Nj=3 QCD 
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Here /,r is the usual parameter of chiral dynam- 
ics, U € SU{3) and 9 is the parameter governing 
CP-violation. At 6 ^ 0,tt there is explicit CP- 
violation whereas at 9 — Sir = tt corresponding 
to negative quark mass spontaneous CP-violation 
is expected. The ^-dependence of the energies of 
different ground states are illustrated by fig. |l|. 

This picture is the solution of the puzzle em- 
phasized in 1^: the vacuum energy in the Ac — > 
oo limit is 
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E{e) = N^h[ — ] =E{0 + 2n) 



The solution can only be a multi-branched func 
tion 



E{e) = min h 

k 



Ac 



(9) 



Numerical simulations could prove this picture 
following from low energy effective chiral La- 
grangians. 
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Figure 1. The energies of different branches of 
ground states as a function of the ^-parameter. 
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